STOCHASTIC LATTICE MODELS WITH SEVERAL ABSORBING STATES 
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We study two models with n equivalent absorbing states that generalize the Domany-Kinzel cellular 
automaton and the contact process. Numerical investigations show that for n = 2 both models be- 
long to the same universality class as branching annihilating walks with an even number of offspring. 
Unlike previously known models, these models have no explicit parity conservation law. 
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I. INTRODUCTION 

The study of stochastic lattice models exhibiting a con- 
tinuous phase transition from an active phase into an 
absorbing state is a field of growing interest. In these 
models the dynamical processes take place close to an 
absorbing state, i.e. a configuration once reached, the 
system cannot escape from. Most of them belong to the 
universality class of directed percolation (DP); the best 
known examples are DP lattice models , the Domany- 
Kinzel cellular automaton the contact process 
Schlogl's first and second model Q , and branching anni- 
hilating walks with an odd number of offspring B . In a 
field-theoretical formulation these models can be related 
to Reggeon field theory |6|, which was proven to be in 
the same universality class as directed percolation 0| . 

The variety of DP models led Janssen and Grassberger 
to the conjecture that in one-component models all con- 
tinuous phase transitions from an active phase to a single 
absorbing state are in the DP universality class [|| . How- 
ever, the known examples for DP include even more com- 
plicated systems, e.g., multi component systems and 
models with several absorbing states Some models 
that were initially thought to be in different classes were 
later found to belong to the DP class as well [0. Thus 
the directed percolation universality class is extremely 
robust and covers a wide range of models. 

Among the models with absorbing states only a few 
exceptions are known that do not belong to the DP uni- 
versality class. During the past few years it became clear 
that they represent a universality class that is different 
from that of directed percolation. The known examples 
are the models A and B of probabilistic cellular automata 
p^ , p^ , nonequilibrium kinetic Ising models with com- 
bined zero- and infinite-temperature dynamics , inter- 
acting monomer-dimer models jl5[] , and branching anni- 
hilating walks (BAW's) with an even number of offspring 
Ip^-^. The common feature of all these models is that 
the number of particles (or kinks) is conserved mod 2. 
Therefore the new class is sometimes referred to as the 
parity- conserving (PC) class. 

Initially parity conservation was thought to be the rea- 
son for the emergence of the different universality class. 
However, Park and Park recently showed that in the 
example of an interacting monomer-dimer model a weak 
parity-conserving external field can force the system back 



to the DP class. They concluded that the essential prop- 
erty of the class is not parity conservation, but a symme- 
try among different absorbing states. 

In order to address this question, we propose two one- 
dimensional models that exhibit a phase transition from 
an active to an inactive phase consisting of n equiv- 
alent absorbing states. These models generalize two 
well-known stochastic models for directed percolation, 
namely, the Domany-Kinzel cellular automaton ||^ and 
the contact process Q . We can conclude from numerical 
simulations that for n = 2 both models belong to the PC 
class. 

The models we define are interesting for various rea- 
sons. As generalizations of well-known DP models they 
give a better physical understanding of models in the 
PC class. Moreover, they are defined by ordinary two- 
site nearest-neighbor interactions (rather than three- or 
four-site interactions) and can be generalized easily to 
both higher dimensions and a higher number of absorb- 
ing states. In addition, unlike previously known models, 
our models do not explicitly conserve parity. This con- 
firms that the symmetry among the absorbing states is 
indeed the only essential property of models in the PC 
class. 

Before defining the models let us present an intuitive 
idea how DP models can be generalized to models with 
several absorbing states. Directed percolation models are 
usually defined on some d-dimensional lattice whose sites 
can be either active (wet) or inactive (dry). If all sites 
are inactive, the system is in an absorbing configuration 
from that it cannot escape. In the presence of active 
sites, the system evolves in time according to specific 
local processes. Although microscopically these processes 
may be defined differently, most of the DP models have 
the feature in common that their time evolution seen on 
a large scale is subject to the following rules: 

(a) Inactive (dry) spots are created randomly within 
active (wet) islands. 

(b) The boundaries between active and inactive 
domains fluctuate in a way that active islands 
are biased to grow. 

Both processes (a) and (b) compete with each other. De- 
pending on their probabilities the system can be in two 
different phases. If the probability for (a) is very small. 
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the system is in the active phase where, starting with a 
nonzero density of active sites, active chisters percolate 
constantly. If the probability for (a) is very large, the 
system is in the inactive phase where active clusters die 
out wherefore eventually the system enters the absorb- 
ing state. At the percolation threshold the system goes 
through a continuous phase transition where the details 
of the local processes become irrelevant and long-range 
fluctuations can be observed. According to the DP con- 
jecture 1^, we assume that all models with a single ab- 
sorbing state defined in the spirit of rules (a) and (b) 
belong to the universality class of directed percolation. 

The main idea of the present work is that a generaliza- 
tion of the above rules to n equivalent absorbing states 
generates universality classes different from DP, in par- 
ticular to the PC class in the case of two symmetric ab- 
sorbing states. Such a generalization can be defined as 
follows. Let us assume that each inactive site carries a 
"color" labeled by 1, . . . ,n. The simplest generalization 
of the rules (a) and (b) is the following: 

(A) Inactive spots of random color 1, . . . , n are created 
randomly within active islands. 

(B) The boundaries between active and inactive 
domains fiuctuate in a way that active islands 
grow. 

(C) Boundaries between inactive domains of different 
colors are not allowed to stick to each other 
irreversibly. They are free to separate again 
leaving active sites in between. 

Rule (A) and (B) are straightforward generalizations of 
(a) and (b). Again both processes compete with each 
other and lead to a phase transition from an active to an 
inactive phase. Rule (C) is different and distinguishes 
the different colors. Roughly speaking, this rule tells us 
that between two inactive domains of different colors a 
thin film of wet sites is preserved. The importance of 
this rule becomes obvious by looking at the contrary: 
If domains of different colors were allowed to stick to 
each other irreversibly, the colors would then be irrele- 
vant. This would mean that the process is compatible 
with the previous rules (a) and (b) and thus belong to 
directed percolation. Rule (C) allows wet sites between 
absorbing domains of different colors to survive for a long 
time. This slows down the relaxation towards one of the 
absorbing states and therefore we expect systems with 
several absorbings states to be "more active" than usual 
DP models. 

Another important requirement is that the rules are 
symmetric under global permutation of the colors. We 
will show that if this symmetry is broken, one of the 
colors begins to play a dominant role so that the phase 
transition is again in the DP universality class. Although 
these rules give only an intuitive description rather that 
a strict definition, they will help us to define models with 
several absorbing states. 



II. DEFINITION OF THE MODELS 

A. Model I: Generalized Domany-Kinzel cellular 
automaton 

In the Domany-Kinzel model the state at a given 
time t is specified by binary variables {si}, which can 
have the values A (active) and / (inactive). At odd (even) 
times, odd- (even-) indexed sites are updated according 
to specific conditional probabilities. This defines a cellu- 
lar automaton with parallel updates (discrete time evo- 
lution) acting on two independent triangular sublattices: 




FIG. 1. Update in the Domany Kinzel model. 



The conditional probabilities in the Domany-Kinzel 
model P{si^t+i \ Si-i,t , Si+i,t) are given by 

P{I\I,I)^l, (1) 
P{A\A,A)^q, (2) 
P{A\I,A)^P{A\A,I)^p, (3) 

and P(/|sj_i, s,j+i) + P{A\si-i, Si+i) = 1, where < 
p < 1 and < q < 1 are two parameters. Equation 
(|l|) ensures that the configuration ...,/,/,/,... is the 
absorbing state. The process in Eq. (|^) corresponds to 
rule (a) and describes the creation of inactive (dry) spots 
within active (wet) islands with probability 1 — q. The 
random walk of boundaries between active and inactive 
domains is realized by the processes in Eq. (|^) . According 
to rule (b), DP transitions can be observed only if p > 
^ when active (wet) islands are biased to grow. The 
processes and their probabilities can be summarized in 
the form of a probability table: 



Si-l, Si+i 


P(A|s,j_i,Si+i) 


P(/ 1 Sj_l, Sj+l) 


AA 


q 


l-q 


AI 


p 


1-p 


lA 


p 


l-p 


II 
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We now define a generalization of the Domany-Kinzel 
model following the rules (A)-(C) (hereafter referred 
to as model I). This model has n + 1 states per 
site: one active state A and n different inactive states 
/i , /2 ■ The conditional probabilities are given by 

(kj = 1, . . . k ^ I) 

P{Ik\h,h) = l. (4) 
P{A\A,A) = l-nP{h\A,A)=q, (5) 

PiA\h,A)=P{A\A,h)=Pk, (6) 
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P{h\h,A) 



P{A\hJi] 



Pk . 



(7) 



where we study the symmetric casepi, . . . ,pn = P- Equa- 
tions (||)~(|^) are straightforward generalizations of Eqs. 

"^^^ '^^^y different process is the creation of ac- 
tive sites between two inactive domains of different colors 
in Eq. (|^) according to rule (C). For simplicity we chose 
the probability of this process to be equal to one. We 
may also use a probability less than one, but it turned 
out that this does not change the critical properties of 
the system. 

For n = 1 the model defined above reduces to the origi- 
nal Domany-Kinzel model. In Sec. |l| we will investigate 
the generalized Domany-Kinzel model with n = 2 absorb- 
ing states. The corresponding probability table reads: 





P{A\SI,S2) 


P{h\si,S2) 


P{h\si,S2) 


AA 
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1-9/2 


1-9/2 


Ah 
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1-p 





Ah 


p 





1-p 


hA 
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1-p 





hA 


p 





1-p 


hh 





1 





hh 


1 








hh 


1 








hh 
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B. Model II: Generalized contact process 

The one-dimensional contact process is the simplest ex- 
ample for a DP model with continuous time evolution . 
Its dynamics is defined by nearest-neighbor processes 
that occur spontaneously due to specific rates (rather 
than probabilities). In numerical simulations models of 
this type are usually realized by random sequential up- 
dates. This means that a pair of sites {si, Si+i} is chosen 
at random and an update is attempted according to spe- 
cific transition rates w{si^t+dt , Si+i^t+dt \ Si^t , Si+i,t). 



— O- 




t-Hlt 



FIG. 2. Update in the contact process. 

Each attempt to update a pair of sites increases the 
time t hy dt = 1/N, where N is the total number of 
sites. One time step (sweep) therefore consists of N such 
attempts. The contact process is defined by the rates 

w{A,I\A,A) ^w{hA\A,A) ^ X, (8) 
wihI\A,I)=w{hI\hA)=^l, (9) 
w{A,A\A,I) ^wiA,A\hA) = 1, (10) 



where A > and fi > are two parameters (all other 
rates are zero). Equation (^ describes the creation of 
inactive (dry) spots within active (wet) islands corre- 
sponding to rule (a). Equations (^) and (^^ describe 
the shrinkage and growth of active islands according to 
rule (b). In order to fix the time scale, we chose the rate 
in Eq. (|lO|) to be equal to one. The active phase is re- 
stricted to the region fi < 1 where wet islands are likely 
to grow. 

As in the case of the Domany-Kinzel model, we define 
a generalization of the contact process by introducing n 
different inactive states h, h, ■ ■ ■ , h- The dynamics of 
the generalized model (model II) is defined by the rates 

wiA, h \A,A)= w{h,A I A, A) - \/n , (11) 

w{h ,h\A, h) ^w{h, h I h , A) ^ , (12) 

w{A,A\A,h) ^w{A,A\h,A) ^1, (13) 

w{h,A I /fe, Ii) = w{A, h I h,Ii) - 1 , (14) 

where fc, Z = 1, . . . , n and k ^ I (all other rates are zero). 
Again we consider the symmetric case /xi, . . . — fi. 
Eqs. (|ll|)-(|l|) are generalizations of Eqs. (|)-(|l^). The 
new rule (C) is implemented by Eq. ( |l4[) which describes 
the creation of active sites between two inactive domains 
of different colors. For n = 1 the model defined above is 
reduced to the usual contact process (P)-(p^). 



Model I 



Model II 
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inactive (dry) 

n=2 



active 
(wet) 



n=l 



inactive (dry) 



active (wet) 



n=2 : 



n=l 



FIG. 3. The phase diagrams of models I and II for two 
absorbing states. The dashed lines indicate the corresponding 
transition lines for directed percolation. The explanation of 
the points A,B, and C can be found in the text. 



C. Phase diagrams 

The phase diagrams of both models are shown in Fig. ||. 
The active (inactive) phase is characterized by a nonzero 
(vanishing) density of active sites in the thermodynamic 
limit. Both phases are separated by a phase transition 
line (the bold line in Fig. ||). The dashed line indicates 
the corresponding phase transition for a single absorbing 
state. Comparing both lines we notice that generally 
models with two absorbing states tend to be more active 
than their DP counterparts (for exceptions see Ref. |^ ) . 

We checked numerically that as in DP the critical ex- 
ponents of the generalized models are the same all along 
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the phase transition hne. The only exceptions are the 
ending points {p,q) — (1/2,1) and (A,/x) = (0,1), where 
the transition Unes for n — 1 and n = 2 intersect (marked 
by A in Fig. ||) . In these points rule (A) is no longer valid 
and the entirely active configuration . . . AAA . . . emerges 
as an additional absorbing state. This leads to a different 
universality class, which, in the case of n = 1, is referred 
to as compact directed percolation. 



still be seen o n larg e scales. We will return to this obser- 
vation in Sec. IV C . 



III. TWO SYMMETRIC ABSORBING STATES: 
NUMERICAL RESULTS 

A. Monte Carlo simulations 



D. General properties 

A typical time evolution of models with two absorbing 
states is shown in Fig. ^. In the active phase (A < Ac) 
small inactive islands of random color are generated and 
exist only briefly. Approaching the phase transition their 
size and lifetime grows while the density of active sites 
decreases. Notice that according to rule (C) a thin film 
of active sites separates different inactive domains. 

An important property of models with several absorb- 
ing states is a very different relaxation towards the ab- 
sorbing configuration. For DP in the inactive phase the 
order parameter p is known to decay exponentially in 
time. However, this is not true for models with two 
absorbing states. As shown in Fig. ^ (A > Ac), start- 
ing from a random initial configuration, large domains of 
different colors are formed. These domain walls survive 
and diffuse until they annihilate mutually. In this an- 
nihilation process the density of domain walls is known 
to decay algebraically like p{x) ~ Because of 

the slow relaxation numerical simulations of models with 
several absorbing states are more difficult to perform. 




X<Ac A.=A,c 



FIG. 4. Simulation of model II for n — 2 starting from a 
random initial condition. The two different types of inactive 
domains are shown in black and grey. The active sites between 
the domains are represented by white pixels. 

It should be emphasized that in the models defined 
above there is no explicit parity conservation on the mi- 
croscopic level: In each local update no more than one 
site is modified [cf. Eqs. ([l|)~(|l|)]. Therefore it is 
impossible to create more than two kinks or particles 
per update (a nontrivial paritiy-conserving dynamics re- 
quires the generation of at least three kinks X 3X, 
2X — > 0). Nevertheless, the annihilating domain walls 
described above, by their very nature, obey a parity- 
conserving dynamics. Therefore parity conservation can 



In order to measure the critical exponents of models 
I and II in the case of two absorbing states, we perform 
dynamic Monte Carlo simulations (see, e.g., fl^). We 
use defect dynamics, i.e., we start with an initial config- 
uration where all sites are in the inactive state Ii except 
for one active site in the center. The system then evolves 
along the dynamic rules of the model. For various values 
of the parameters near point B in Fig. || we perform 10^ 
independent runs up to 5000 time steps. However, most 
of them stop earlier because the system enters into the 
absorbing configuration where all sites are in the state 
In order to avoid finite-size effects, we adjust the sys- 
tem size after each time step according to the actual size 
of the active cluster. As usual in this type of simulations, 
we measure the survival probability P{t), the number of 
active sites N(t), and the mean square of spreading from 
the origin R^(t) averaged over active runs. At criticality, 
these quantities are expected to scale algebraically in the 
long-time limit 

p(^)-^-^ Ar(t)-i\ R^{t)r^e'. (15) 




0.0005 0.0010 0.0015 0.0020 J/^ 0.01 0.015 0.02 0.03 



FIG. 5. Numerical results for model I. The effective expo- 
nents 5{t), Tj{t) and z{t) are obtained from dynamic simula- 
tions for p = g = 0.5665, 0.5670, . . . , 0.5685. The densitiy of 
active sites p(e) is measured in static simulations for different 
values of e — p — pc- The slope of line in the log-log plot gives 
an estimate for the exponent /3. 



*In an infinite system, there is no way to reach the other 
absorbing configuration /2. 
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0.0005 0.0010 0.0015 0.0020 -j^^^ 0.01 0.015 0.02 0.03 0.05 0.07 g 



FIG. 6. Results of analog simulations of model II. The pa- 
rameters vary in the range A = /i = 0.624, 0.626, . . . , 0.632. 

The critical exponents are related to the exponents /3, 
v^, and i^ii by 



<5 = 



2v^ 



(16) 



Thus, from our numerical results we can conclude that 
for n = 2 both models belong to the PC universality 
class. 



Known models: 




5 


V 


z 




A and B p 


L 


0.27(8) 






0.6(2) 


BAW n = 2 


121 


0.283(16) 


0.272(12)* 


1.11(2) 


0.94(6) 


BAW n = 2 


IE] 


0.285(2) 


0.000(1) 


1.141(2) 


0.92(3) 


BAW 71 = 4 


IS 


0.286(2) 


0.000(1) 


1.147(4) 


0.922(5) 


kinetic Ising 


14] 


0.27(2) 


0.30(2)* 


1.14(2) 


0.80(8) 


dyn. BAW |22(1 


0.287(1) 


0.000(3) 


1.155(5) 




MDM [151 


0.29(2) 


0.00(2) 


1.34(20) 


0.88(3) 


Present work: 


model I 


0.285(10) 


0.00(1) 


1.15(1) 


0.90(5) 


model II 


0.29(1) 


0.00(1) 


1.15(1) 


0.93(5) 



TABLE I. Critical exponents of models in the PC class. The 
star indicates values measured in kink dynamics where a 
different scaling relation holds. 



B. Symmetry breaking field 



and obey the scaling relation 

4:6 + 27] = dz. 



(17) 



The quantities (jl^) show straight lines in double logarith- 
mic plots. Off criticality, the lines are curved. In order 
to get precise estimates for the scaling exponents, it has 
been useful to consider the local slopes of the curves by 
introducing the effective exponents 



Sit) 



logio(P(i)/P(t/6)) 



login b 



(18) 



and similarly ri{t) and z(t), where logj^g ^ the distance 
used for estimating the slope. Choosing 6 = 5, we mea- 
sured the effective exponents of both models for various 
values oi q = p and X = v. The results of our simu- 
lations are shown in Figs. || and ^. Off criticality, the 
curves for S(t) and r]{t) show negative or positive cur- 
vature. The figures give us an estimate of the critical 
points Pc = 0.5673(5) for model I and Ac = 0.628(2) for 
model II. The estimates for the critical exponents are 
S = 0.285(10), Tj = 0.00(1), and z = 1.15(1) for model 
I and S = 0.29(1), = 0.00(1), and z = 1.15(1) for 
model II. 

The exponent (3 has been obtained directly in static 
simulations by measuring the steady state density p in 
the active phase near the critical point. Although this 
method is known to be quite inaccurate to determine (3 
and we measured only over one decade in e, we found the 
reasonable values jS = 0.90(5) for model I and (3 — 0.93(5) 
for model II. 

The estimates of the critical exponents agree with pre- 
vious results for models in the PC class (cf. Table I). 



Recently Park and Park showed in the example of an 
interacting monomer-dimer model that if the symmetry 
among the absorbing statess is broken by an external 
field, the DP universality class is recovered In or- 
der to verify this observation, we introduce an external 
field by modifying the growth rates for inactive islands of 
different colors. This can be done by choosing different 
Pk in Eq. (^) [fik in Eq. (p^)]. Because of the different 
growth rates, one of the colors is going to play a dom- 
inant role wherefore in the large-scale limit the system 
evolves as if it had only a single absorbing state. 




0.02 0.04 0.06 0.08 \/l 



0.01 0.015 0.02 0.03 £ 



FIG. 7. Symmetry breaking field: Simulation of model I for 
slightly different growth rates (pi,2 ~ g±0.02). DP exponents 
are recovered. 
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To demonstrate this, we repeat the above simulations 
of model I forpi_2 = g±0.02. The critical point is shifted 
to qc — 0.6245(10). The exponents we measure (cf. Fig. 
I) arc S ^ 0.16(2), 77 = 0.32(4), z = 1.27(4), and /3 = 
0.28(1), which agree with the DP exponents S ~ 
0.159(1), 77 = 0.314(1), z = 1.266(2), and /? = 0.2765(1). 

If the symmetry of a model with n > 2 absorbing states 
is partially broken, a subset of m colors starts to play a 
dominant role. We expect that such a system behaves at 
criticality like a model with m absorbing states. 

IV. OTHER SPECIAL CASES 
A. Compact clusters 

As mentioned before, the Domany-Kinzel model and 
the contact process have a line in their phase diagram 
where the entirely active configuration emerges as an 
additional absorbing state. On this line active clusters 
are not fractal but compact, wherefore it is called the 
compact directed percolation line. Here the dynami- 
cal processes are exactly solvable by reducing it to an 
annihilation-diffusion process of kinks 2X 0. Al- 
though this leads to a different universality class at the 
phase transition point (point A in Fig. |^), compact DP 
has been used in many cases to improve the understand- 
ing of ordinary DP. 

A similar situation exists in generalized models with 
several absorbing states. The fact that the transition 
A point is identical to that of ordinary compact DP in- 
dicates that this point is exactly solvable in all cases. 
However, the translation into a kink language is slightly 
more complicated. Since kinks between inactive domains 
of different colors cannot exist [rule (C)], only n types of 
kinks Xk between active and inactive domains {Alk, Ik A) 
play a role. These kinks have to occur in pairs and un- 
dergo an annihilation-diffusion process 2Xk — > (kinks 
of different types cannot annihilate). It is important that 
there is no more generation of randomly colored inactive 
domains. Therefore dynamical simulations such as those 
in the previous section starting from an initial configura- 
tion with only one type of inactive sites yield the same 
results for all n, i.e., (5=1/2 and z = 1. However, com- 
pact percolation processes are known to depend strongly 
on initial conditions so that generally the situation may 
be more complicated. 

B. Simulations in higher dimensions 

The models presented in this paper can easily be gen- 
eralized to higher dimensions. This is particularly simple 
in the case of model II since its definition (|ll|)-(|ll) 

can 

be used in any dimension. In simulations we observed 
that for d = 2 this system has a phase transition al- 
though the relaxation towards one of the absorbing states 



in the inactive phase is extremely slow. Figure g shows 
typical configurations in the active phase near critical- 
ity. While in ordinary DP active clusters are separated 
spatially, active sites in two-dimensional models with two 
absorbing states are arranged in fractal "lines" along the 
boundaries of inactive islands. Repeating the simula- 
tions described above on a 80 x 80 lattice we obtained 
Ac = 0.99(1), S = 0.9(1), 7j = 0.00(5), and z = 1.0(1). 
These results agree roughly with the mean-field expo- 
nents (5 = 1, 7y = 0, z = 1 and f3 = 1 (see next Section). 
Therefore, we conjecture that 1 < c?c < 2 is the upper 
critical dimension of systems with two absorbing states. 




n=l (DP) n=2 (PC) 



FIG. 8. Simulation snapshot of two-dimensional systems 
with one and two absorbing states in the active phase near 
criticality. White dots denote active sites. 

C. Relation to BAW's 

We already mentioned that models I and II have the 
same critical behavior all along the phase transition line 
(except point A in Fig. ||). Moving along this line we 
can control the mean size of active islands, which is in- 
finite at point A, a few sites at point B, and one site at 
point C (in model II the latter case corresponds to tak- 
ing A 00). In one dimension point C can be related to 
branching walk models since each active site can be inter- 
preted as a single walker. For n = I these walkers diffuse 
and interact by 2X ^ X , X ^ 0, which results in a DP 
process. For n > 2 the situation is more complicated. 
Let us denote by Xj^. a walker separating inactive do- 
mains Ij and Ik- Then the processes are Xjk <-> XjiXik 
(with random I) and Xjj — > 0. These processes clearly 
do not conserve parity. However, because of Xjj — > 
only walkers between domains of different color j ^ k 
survive for a long time [cf. rule (C)]. Hence for n = 2 
the majority of walkers react like 

X12 AriiAri2 X12 

X12 AriiAri2 ^i2-''^2i-''^i2 

with analog reactions for A'21 

Thus, in a long time limit the walkers undergo an effective 
reaction of the type X 3X, 2X 0, which is a BAW 
with two offspring in one dimension. 
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The relation to BAW modes is even more general and 
holds not only at point C but everywhere on the phase 
transition line except for point A. The walkers then have 
to be identified with domain walls seperating inactive 
domains of different colors. At first glance this seems to 
be contradicting: Domain walls, by their very definition, 
obey a local parity-conserving dynamics. On the other 
hand, it is obvious that the local processes do not con- 
serve parity. However, domain walls in our models are 
extended objects. Their thickness fluctuates and varies 
along the phase transition line from typically one site at 
point C up to infinitely many sites at point A. It is im- 
portant to notice that the domain walls simply cannot be 
identified with active islands (which may also occur be- 
tween inactive domains of the same color) but require a 
more complicated definition. Although there is no micro- 
scopic parity conservation, a careful analysis shows that 
the dynamical rules ensure that all microscopic processes 
violating parity conservation have a very short lifetime. 
This is the reason why an effective parity -conserved dy- 
namics can be recovered in the limit of large scales in 
time and space. 

In higher dimensions d > 2 the physical properties of 
BAW's with an even number of offspring are governed 
by the mutual annihilation of the walkers However, 
the models presented in this paper behave very differ- 
ently in higher dimensions, which makes it impossible to 
relate them to BAW's. As shown in Fig. ^ the active 
sites in models I and II arrange themselves as {d — 1)- 
dimensional surfaces separating inactive domains of dif- 
ferent colors. Thus they cannot be interpreted as point- 
like random walkers. In d > 2 dimensions we therefore 
expect BAW's with an even number of offspring to be in 
a different universality class from the present models. 



D. More than two absorbing states 

No results were obtained for n > 3 symmetric absorb- 
ing states in one dimension. It turned out that it is im- 
possible to determine the critical point because the plots 
for S{t), rj{t), and z[t) show only one type of curvature. 
This observation agrees with recent results obtained for 
an iV-species generalization of BAW models with an even 
number of offspring . For N >2 these models are al- 
ways in the active phase and their critical exponents are 
described by yet another universality class. 



E. Mean-field approximation 

Denoting by pk the density of inactive sites Ik, the 
processes in model II imply the mean-field equation 



d_ 

di 



2A 



Pk 



pI + PPoPk + pI- Pk 



(19) 



means that in mean-field models with more than two 
absorbing states are always in the active phase (see also 
Ref. 1^). For all n the density of active sites po scales 
like (A — Ac)'' with the mean-field exponent (3—1. This 
is close to the measured values (3 ^ 0.9 in the model with 
two absorbing states, which indicates that we are already 
close to the upper critical dimension. 



V. SUMMARY AND DISCUSSION 

We have shown by the example of the Domany-Kinzel 
model and the contact process that lattice models for di- 
rected percolation can be generalized to models with n 
symmetric absorbing states. Numerical simulations lead 
to the conclusion that such models in one dimension with 
two absorbing states belong to the PC universality class. 
Since these models do not explicitly conserve any quan- 
tity mod 2 they show that rather than parity conservation 
the symmetry among the absorbing states is the origin 
for the emergence of a different class. As soon as this 
symmetry is broken, the critical behavior jumps back to 
DP. 

The symmetry used in our models is the group of per- 
mutations Sn- No reliable numerical results could be ob- 
tained for S3 and higher symmetries. However, one may 
also introduce other symmetries such as cyclic groups 
(e.g., Z3) and investigate whether they define different 
universality classes (cyclic symmetries of this type ap- 
pear, e.g., in the three-candidate voter model p5| ). It 
would be also interesting to examine models with more 
than one symmetric active states. 

For a better understanding it would be desirable to 
find an appropriate field-theoretical description of the 
model. For BAW's with an even number of offspring 
this has been done recently in Ref. There arc sur- 

prising results, in particular one has two different critical 
dimensions, one of them (dci = 4/3) related to the prop- 
erties of active clusters and the other (dcj = 2) related 
to the annihilation process. Although the application 
of this theory to the present type of model may not be 
transparent, Cardy and Tauber were able to identify an 
S2 symmetry on an operator level. This again indicates 
that this symmetry plays an important role in the PC 
class. 

A field theory for the present type of model should 
be different from that of BAW's because of its different 
phen omenology in higher dimensions (see Sees. IV B| and 
IV C). Nevertheless, both theories should give identical 
results in one dimension. The simplest ansatz for such a 
field theory would be to add a diffusion like term and n 
noise fields to the mean-field equations 



where pQ — 1 ~ X]J=iPj- Choosing A = /i, the critical 
point is Ac = 1 for n — I and Ac = 00 for 71 > 2, which 



d 2A 

— Pfc(a;, t) = — pI{x, t) + ppo{x, t)pk{x, t) + pl{x, t) 
at n 

-Pk{x, t) + D V^pk{x, t) + 77fc(a;, t) . (20) 
However, even if this were correct, the derivation of the 
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correlations in the noise remains a highly nontrivial prob- 
lem. Therefore, the development of appropriate field 
theories is a challenge towards a better understanding 
of universality classes appearing in systems with several 
absorbing states. 
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